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Abstract 

Wc introduce a now micro/macro decomposition of coUisional kinetic equations which naturally incorporates the 
exact space boundary conditions. The idea is to write the distribution fonction / in all its domain as the sum of a 
Maxwellian adapted to the boundary (which is not the usual Maxwelliaii associated with /) and a reminder kinetic 
part. This Maxwellian is defined such that its 'incoming' velocity moments coincide with the 'incoming' velocity 
moments of the distribution function. Important consequences of this strategy are the following, i) No artificial 
boundary condition is needed in the micro/macro models and the exact boundary condition on / is naturally 
transposed to the macro part of the model, ii) It provides a new class of the so-called 'Asymptotic preserving' 
(AP) numerical schemes: such schemes are consistent with the original kinetic equation for all fixed positive 
value of the Knudsen number e, and if e — >^ with fixed numerical parameters then these schemes degenerate 
into consistent numerical schemes for the various corresponding asymptotic fluid or diffusive models. Here, the 
strategy provides AP schemes not only inside the physical domain but also in the space boundary layers. We 
provide a numerical test in the case of a diffusion limit of the one-group transport equation, and show that our 
AP scheme recovers the boundary layer and a good approximation of the theoretical boundary value, which is 
usually computed from to the so-called Chandrasekhar function. 

Resume 

Nous introduisons une nouvelle decomposition micro/macro pour les equations cinetiques coUisionnelles qui tient 
compte de fagon exacte des conditions aux bords en espace. L'idee est de decomposer la fonction de distribution / 

en une Maxwcllicnnc adaptee au bord (qui n'est pas la Maxwellienne usuellc associoe a /) et une partie cinetique 
restante. Cette nouvelle Maxwellienne est definic dc sorte que ses moments 'entrants' en vitcsse soicnt les memes 
que les moments 'entrants' de /. Des consequences importantes de ccttc strategic sont les suivautos. i) Aucune 
condition artificielle n'est necessaire et les conditions de bord sur / sont naturellement transposees sur la partie 
macroscopique de la decomposition, ii) La strategic produit une nouvelle classe de schemas dits 'Asymptotic 
preserving' (AP) : ccs schemas sont consistants avec le modcle cinetique original pour toutc valeur positive du 
nombre de Knudsen £, et si e — >■ et les parametres numeriques sont fixes, alors ces schemas degenerent en des 
schemas qui sont consistants avec les divers modeles asymptotiques associes. Ici, les schemas obtenus reproduisent 
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tres bien les couches limites diffusives. Nous donnons un resultat numcriquc dans le cas d'une limite de diffusion 
et montrons que la valeur au bord theorique (calculee habituellement a partir de la fonction de Chandrasekhar) 
est bien approchee numeriquement. 



Version frangaise abregee 

Le but dc cc travail est d'introduire une nouvellc decomposition micro/macro pour les equations 
cinetiques, capable par construction d'integrer les conditions aux bords exactes du modele. Dans les 
approches usuelles, les modele macrosopiques (asymptotiques ou aux moments) utilisent generalement 
dcs moyenncs en vitcssc de la fonction / ainsi que leurs flux, et les valours des ces moyennes au bord 
ne peuvent etre en general deduites de la condition au bord sur la fonction distribution /. En effet, si 
la condition au bord en espace est de type Dirichlet sur /, alors la valeur de / au bord ne peut etre 
imposce que pour les vitesses entrantes. Par consequent les moyennes (sur tout le domaine des vitesses) 
de / ne sont pas donnees a priori, et des calculs supplementaires (probleme de Milne) et/ou des conditions 
artificielles sont en general necessaires pour bien poser le modele macroscopique. 

La strategic presentee ici permet dc rcsoudre ce probleme et constitue une base pour le dcvcloppement 
des schemas dit 'Asymptotic Preserving' (AP) qui ont en plus la specifite importante de pouvoir integrer 
les conditions aux bords de fagon exacte. Avant de detailler la strategie dans les sections suivantes, nous 
rappelons brievement la problematique des schemas multi-echelles dans ce contexte. Quand le nombrc 
de Knudsen e devient petit, une raideur apparait dans ces modeles qui fait en sorte que leurs simula- 
tions numeriques par des schemas explicites en temps deviennent rapidement inaccessibles (contraintes de 
type At = 0(e) ou Le probleme est alors de construire des schemas numeriques pour I'equation 

cinetique consideree qui s'affranchissenent de ces contraintes et degenerent en des discretisations consis- 
tantes avec les modeles asymptotiques quand £ ^ 0. Plusieurs travaux ont ete effectues a ce propos, 
voir par exemple [3,4,5] pour les limites de diffusion. Une bibliographic plus complete sera donnee dans 
une version detaillee de ce travail [8]. Une approche generale a ete introduite dans [1,6,7], permettant de 
construire des schemas numeriques qui sont consistants avec les modeles cinetiques etudies pour toutes 
les valeurs de e > 0, et degenerent quand £ — en des discretisations consistantes avec les systemes 
d'Euler, de Navier-Stokes compressibles, ou avec la limite de diffusion. Elle est basee sur la decomposition 
micro-macro qui permet d'ecrire de manierc cquivalente I'equation cinetique de depart sous la forme d'un 
systeme couplant une equation fluide ou dc diffusion avec une equation sur la partie cinetique restante. 

Le but ici est de completer ccttc scrie de travaux sur les schemas micro/macro, en developpant une 
nouvelle strategie qui tient compte des conditions de bord en espace de fagon exacte. Elle consiste a choisir 
une partie macrosopique 'integrce sur les vitesses entrantes', qui est compatible avec les conditions aux 
bords, et de la coupler avec la partie cinetique restante. En plus d'integrer naturellcmcnt les conditions 
aux bords, cette strategie constitue une base pour la construction des schemas numeriques qui ont toutes 
les bonnes proprictcs dcs schemas micro/macro decrits ci-dessus. Cette note donne les grandes ligncs dc 
la strategie qui sera dcveloppee en detail dans [8]. Nous donnons aussi les resultats d'un premier test 
numerique dans le cas d'une limite diffusive, qui montre que le schema construit est bien AP et permet 
de retrouver la bonne couche limite au bord. Des experiences numeriques plus completes seront effectuees 
dans [8]. 
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1. Boundary matching micro/macro formulations of kinetic equations: the general problem 

It is a usual challenge to design efficient numerical schemes for kinetic equations which are consistent 
with the kinetic model for all positive value of the Knudsen number £, and degenerate into consistent 
schemes with the asymptotic models (compressible Eulcr and Navicr-Stokcs, diffusion, etc) when e — > 0. 
On this subject, many works can be quoted, see for instance [3,4,5] and a more complete bibliography in 
the forthcoming paper [8]. Recently [1,6,7], a general strategy based on the micro/macro decomposition 
was introduced. It consists in writing the kinetic equation as a system coupling a macroscopic part (say 
a Maxwellian) and a reminder kinetic part. The main advantage of this method is its robustness and 
easy adaptability. However, the general problem of transfering boundary conditions from the distribution 
function / to the two parts of the micro/macro systems has not been solved and artificial boundary 
conditions were needed. 

Our aim is to develop a new micro/macro decomposition of collisional kinetic equations which naturally 

incorporates the exact space boundary conditions (BC). In fact, it is well known that during the derivation 
of a macroscopic model from kinetic equation (via moments or Chapman- Enskog approach), the exact 
boundary conditions on / are generally lost. Indeed, if one takes moments of the kinetic equation, then 
the values of the fluxes at the boundary are required. But, in general (for inflow BC for instance), 
the distribution function / is known on the boundary for only incoming velocities and, therefore, their 
moments are not prescribed at the boundary. 

The idea is to decompose the distribiition fonction / in its domain as the sum of a Maxwellian part 
adapted to the boundary (which is not the usual Maxwellian associated with /) and a reminder kinetic 
part. This Maxwellian is defined on the whole domain such that its 'incoming' velocity moments coincide 
with the 'incoming' velocity moments of the distribution ftmction. Important consequences of this strategy 
are the following, i) No artificial boundary condition is needed in the micro/macro models and the exact 
boundary conditions on / are naturally shared by the macro part and the kinetic part. Note that the 
traditional approaches cannot provide exact boundary conditions for the asymptotic and kinetic models 
separately, and extra calculations (Milne problems for instance) or artificial boundary conditions are 
necessary in general, ii) It provides a new class of the so-called 'Asymptotic preserving' (AP) numerical 
schemes: such schemes are consistent with the original kinetic equation for all fixed positive value of the 
Knudsen number £, and when £ — )■ with fixed numerical parameters, these schemes degenerate into 
consistent numerical schemes for the various corresponding asymptotic models. A further fundamental 
property of the strategy developed here is that it provides AP schemes not only inside the physical domain 
but also in the space boundary layers with exact boundary conditions. We shall refer to this property as 
'Boundary Asymptotic Preserving' (BAP). At the end of this Note, we provide a numerical test which 
illustrates this property. We emphasize that our scheme does not use any artificial boimdary condition: 
it provides at the diffusion limit a good approximation of the theoretical boundary value derived from 
the Chandrasekhar function, without injecting this value in our scheme. This Note summarizes the main 
lines of the strategy and announces a more complete presentation [8], where numerical discretizations of 
this strategy will be constructed and implemented with a large variety of numerical tests. 

2. The case of the diffusion scaling 

Let be a domain in the position (physical) space K'^, with boundary dO., and V he a. domain in the 
velocity space K'', endowed with a measure dfi. We consider a function ui : ft x V ^ R such that for all 
X e dSl and v gV, u){x, v) has the same sign as n{x) ■ v, where n{x) is the outgoing normal vector to dCl 
at X. Let us give explicit examples of oj{x, v) for specific geometries. For a ball centered at the origin, one 
can take ui{x,v) = x ■ v. For a half plane xi > 0, one can choose U!{x,v) = (— i;i,0, ...,0). In dimension 
one, for the interval (0, 1), one can take uj{x,v) = {2x — l)v. We then consider the following transport 
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equation in a diffusive scaling 

edtf + v-S/^f=hf, t>0, {x,v)enxV, f\t=o = finit (1) 

where / is the distribution function of the particles that depends on time t > 0, on position x G il, and on 
velocity v G V. The linear operator L acts on the velocity dependence of / and describes the interactions 
of particles with the medium. Wc assume that there exists a positive equilibrium function £ = £{v) with 
{£) := Jy £dfi = 1 and satisfying {v£) = 0, and that the collision operator L is non-positive and self- 
adjoint in L'^{V,£^^ dp); with null space and image given by Af{L) = Spanjf } = {/ = pf, where p := 
(/) := Jy fdp}, n{L) = {M{L))^ = {/ such that (/) = 0}. When e goes to in (1), it is easily seen 
that / converges to an equilibrium state /o = po{t, x)£{v). The diffusion limit is the equation satisfied by 
the density po and is classically given by 

^tpo + • (kV^Po) = 0, with K={vL-\v£)). (2) 

Following [6] and [7], equation (1) can be written in a micro-macro equivalent form via the decom- 
position f = p£ + g, with p(t,x) = (/) and g = f — p£ {g is not necessarily small as in [6]). Let 11 
be the orthogonal projector in L^{£~^dp) onto the nuUspace of L: = {(j))£. Then, inserting this 
decomposition into the kinetic equation and applying 11 and 7 — 11 successively, one gets after direct 
computations 

dtp+^V,-{vg) = Q, dtg+^^{I-Il){vVxg) = ^[Lg-e£vV^p)\. (3) 

We now emphasize that in this formulation, the space boundary condition on p and g are not known 
because they cannot be inferred from the boundary conditions on / in general. Indeed, in the typical case 
of incoming boundary conditions, we impose 

/(i, X, v) = fb{t, X, v), Vi > 0, V(a;, v) edflxV such that ij{x, v) < 0. (4) 

ft is therefore clear that the values of p{t, x) = Jy fit, x, v)dp cannot bo determined a priori on the 
boundary dQ, since / is only known for incoming velocities v, i.e. such that co{x,v) < 0. Our aim in 
this work is to develop a new micro/macro decomposition which is able to incorporate the boundary 
conditions in an exact way. To this purpose, let us introduce some few further notations: 

V_{x) = {vgV, io{x,v) < 0}, V+{x) = V\V_{x). (5) 

The idea is now the following. Instead of looking for an equation on p as usual, we seek an equation on 
the following 'boundary matching' density and perform the corresponding micro/macro decomposition: 

{f{t,X,-))y_ f 

^(*'^) = 7fv7 ^\ — ' (-)v-= • t'M, f = P^ + 9- (6) 

{C(T,X,.))y_ Jy_ 

When e ^ 0, we know that the solution / of (1) is (at least formally) close to p£ except in initial or 
boundary layers. Therefore the 'boundary matching' density p will be close to p for small s. This shows 
that (6) still a decomposition of / into an asymptotic part (macro part) and a kinetic part (micro part). 
In order to derive the system of equations satisfied by p and g from (1), we first integrate (1) on V- and 
get 

dtp+-{v£)y ■V^p+-{v-V^g)y ^-{Lg)y . (7) 
Substracting this from the equation on /, we obtain the equation on g: 

dtg +-^(y-Vxg-{v- V,g)^_) + J (^^ " ' = ^ {^9 - • (8) 

System (7)-(8) can be replaced by a the more convenient (and still equivalent) system in terms of p = 
(/) = Jy fdp and g = f — p£ as follows: 
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dtp + - {v- V^g) = 0, 

dtg + - • '^^ff - (^^ • '^xg}v_) + ~[^^ ~ ("^)\/_ j • "^^P = ^ (-^5 - {Lg)y j . 

Note that p is linked to p and g by the relations p = p — (g) ~ p — {9)v+ ^^'^ f ~ P ^ is)v+ ^~ 

One main interest of this new miero/macro formulation (9) for the original kinetic equation (1) is 
the following: the fluxes involved in (9) only concern the quantities g or p, and not p. This means that 
numerical schemes of this formulation would only need the values of p and g at the space boundary which 
are completely known from the original boundary condition (4) on /. 



3. The case of a fluid scaling 

We now consider non linear kinetic equations with a fluid scaling 

dtf + v-WJ=-Qif,f), t>0,{x,v)GnxV, f{t = 0,x,v) = finitix,v), (10) 

where C M'^ and V = R'^. The collision operator Q is a quadratic operator acting only on the velocity 
dependence of the distribution function /. Fundamental examples are the well known Boltzmann kernel 
for rarefied gases and the Fokker-Planck-Landau operators for plasmas. In all what follows, we use the 

notations 

miv) = il,v,\v\y2f, and (g) ^ f g{v) dv (11) 

for any scalar or vector function g = g{v). It is well known that the Boltzmann and Landau opera- 
tors Q{f- /) have important physical properties of conservation and entropy: for all / > 0, we have 
{mQ{f, /)) = and {Q{f, /)log(/)) < 0. It is also well known that the equilibrium functions (/ such that 
Q{f, f) — 0) are MaxwcUians: 

M{U){v) = p{2iTT)-^''^cxp(^ -\v- u| V2t) , U ^ {p, pu, p\uf jl + {d/'i)pT) ^ {mM{U)) . (12) 

To any distribution function /, we shall associate its Maxwellian M = M[f] — M{U), that is the fonction 
of the form (12) which has the same moments as /: U = (mf) = {mM[f]) = {mM{U)). 

When £ goes to 0, / approaches a Maxwellian. At the first order in £, the solution / approaches that 
of the compressible Euler system of gas dynamics. This is formally obtained by integrating the kinetic 
equation (10) against m(w) and replacing / by its first order approximation: the Maxwellian which has 
the same first moments as /. At the second order in £, a standard Chapman-Enskog expansion leads to 
the compressible Navicir-Stokcs system. 

In order to design numerical schemes which arc able to reproduce these asymptotics, a general class 
of numerical schemes has been introduced in [1] and [7] on the basis of micro/macro decompositions. 
The strategy in [1,7] consists in decomposing f as f — M{U) + g and transform (10) into an equivalent 
system of equations on M and g. Then, discretizations of this formulations lead to a class of AP numerical 
schemes which are consistent with the fluid limit in the asymptotics £ — >■ 0. However, similarly as above 
for the diffusion limit, the exact space boundary conditions are lost in this decomposition, and one 
should compute approximate EC for the asymptotic models. To solve this problem, we introduce here 
a more suited decomposition matching the boundary. We flrst deflne the new moments U by U{t,x) = 
{m{v)f{t,x,v))y where VL is defined in (5), and the corresponding equilibrium M is the (uniquely 

defined, see [8] for details) Maxwellian such that Jy m{v)Mdv = U{t,x). In some sense, this strategy 
has some similarities with the half-moment method developed in [2]. Note that M is of the form (12), i.e. 
M = M{U), but in general U is different from U. In fact, since the domain of integration is now V- and 
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not M'', the relation between the parameters of M and U is not explicit in general. We now introduce the 
following micro/macro decomposition: 

f = M + g, / m{v)Mdv = U{t,x) = m{v)f{t,x,v)dv. (13) 
Jv- JV- 

When £ — )• 0, it is known (at least formally) that the solution / of (10) approaches its classical Maxwellian 
M: f-M = 0{e). This means that Jy_ m{v){f-M) = 0{s), and then from (13), we have Jy_ m{v)(M- 

M) = 0{e). This gives d + 2 relations bcitwcxni the two MaxwcUians M and M, each of these Mawellians 
being completely determined by d + 2 parameters (from (12)). Therefore M — M = 0(e). In some sense, 
this justifies the name 'micro/macro' for this new decomposition. 

In order to transform (10) into a coupled system on M and g, we first introduce a definition. For each 

Maxwellian M, we define 11 as the orthogonal projection on SpanjAf, vM, \v\'^M} in (m 'v_).This 

projector has an explicit expression, see [8]. We now insert the decomposition (13) into (10): 

dtM + v V^M + dtg + vW^g=^^ {L^g + Q{g, g)) , Lj^g = 2Q(M, g). 

Now instead of writing a system on M and g, it turns out to be more convenient to rather write a system 
on the full moments U ~ {"nif) and g. To obtain it, we simply integrate (10) against m{v) on the whole 
velocity domain V = R'^ for the equation on U, and apply / — 11 to (10) to get the equation on g: 

dfU + Va; • ( / vm{v)Mdv\ + V^; • ( / vm{v)gdv\ = 0, 

\Jv.-i J \Jm.a J (14) 

dtg + (7 - n) (t; • V,M) + (7 - n) (^; • V,^) = \{l i^Md + Q{9, 9)) ■ 

Note that, since f = M + g, the parameters {p,u,T) of the Maxwellian M (according to definition (12)) 
are linked to U and g by the simple relation: {'p,'pu,p\u\'^/2 + {d/2)'pT) = U — {m{v)g) . 

On the other hand, this relation relation cannot determine M at the boundary because U is not known 
at the boundary. However, by construction, M is completely computed at the boundary from the relation 
(m(u)M)^ = {m{v)fb)y_ , where is the incoming boundary condition (4). 



4. A numerical test 

In this Note, we propose to validate the strategy in the case of the diffusion limit only. A more detailed 
numerical validation will be done in [8]. We consider the simple framework of one-dimensional spaces 
a; G O = [0, 1] and v £ V = [—1, 1] with dp = ^dv, and we take Lf = \ j^i fdv — f ~ P — f ■ The 
initial data is finit{x,v) ~ and the boundary data are f{t,0,v) = v for v > and f{t,l,v) = 
for t; < 0. On Figure 1, we plot the densities p obtained by the AP micro-macro scheme matching the 
boundary, for the values e e {1,0.5,0.2,0.05,0.0001} and compare them with the reference solutions. 
For e £ {1, 0.5, 0.2, 0.05}, the reference solutions are computed with a highly resolved explicit scheme for 
(1). For s = 0.0001, the reference solution is computed by the diffusion equation (2) with the theoretical 
value of the Dirichlet boundary condition, which is p(0) = 0.71043... at x = 0, calculated thanks to the 
usual Chandrasekhar function, and p(l) = at x = 1. The numerical results show that our scheme is 
AP in the diffusion regime and is able to reproduce efficiently the boundary layer. The present numerical 
experiment confirmes that the scheme is 'BAP' in the sense described above: without introducing any 
artificial boundary condition, we recover both the diffusive regime and the boundary layer. 
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Figure 1. Left: numerical solutions in the kinetic and transition regimes, £ = !,£ = 0.5, £ = 0.2. The density for the 

micro/macro scheme coincides with the reference calculated with a highly resolved explicit scheme. Right: numerical solutions 
in the transition and diffusion regimes. The density for the micro/macro scheme coincides with the reference. For £ = 0.05, the 
reference is calculated with the highly resolved explicit scheme. For e = 0.0001, the reference is calculated by the diffusion 
equation with the Chandrasekhar value at the boundary. Both figures are plotted at time t = 0.4 and all micro/macro 
schemes are used with a uniform grid of 100 points in x. 
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